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EXECUTIVE  SUMMARY 

'  THE  MAGNETIC  FIELDS  OF  A  HORIZONTAL 

ELECTRIC  DIPOLE  IN  A  SEMI-INFINITE  MEDIUM 

OBJECTIVE 

The  objective  of  this  report  is  to  provide  a  means  of  predicting, 
interactively,  the  magnetic  fields  in  a  quasi-static  range  due  to  a 
submerged  horizontal  electric  dipole. 

APPROACH 

Simple  engineering  expression  were  used,  or  developed,  to  predict  the 
magnetic  fields.  These  expressions  were  programmed  in  ANSI  77  Fortran  to 
produce  graphics  displays  on  a  Tektronix  terminal. 

RESULTS 

The  computer  programs  can  be  used  to  find  the  magnetic  fields  along 
certain  chosen  paths  in  the  upper  half  or  in  the  lower  half  medium.  The 
corresponding  magnetic  curves  for  each  path  can  then  be  plotted.  The 
systems  which  predict  the  results  for  DC  sources  have  been  utilized 
extensively  by  the  David  W.  Taylor  Naval  Ship  Research  and  Development 
Center,  Annapolis,  Maryland. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Computer  simulations  are  viable  tools  to  evaluate  magnetic  fields  in  an 
interactive  environment.  In  order  to  fully  analyze  experimental  data,  it  is 
recommended  that  the  models  be  modified  to  include  multi-layer  effects. 

CAUTION:  This  document  was  prepared  as  part  of  a  larger  effort.  The 

contents  should  not  be  taken  out  of  context  of  that  larger  effort. 

vi 


ABSTRACT 


Various  formulae  for  magnetic  fields  of  horizontal  electric  dipoles 
in  a  semi -infinite  medium  have  been  derived,  in  the  interest  or 
providing  the  ability  to  evaluate  the  magnetic  field  strength  interactively 
at  on-site  field  locations,  reduced  expressions,  which  are  valid  in  the 
quasi-static  range,  for  the  magnetic  field  were  used  in  developing 
computer  generated  plots. 

Existing  expressions  were  used  to  plot  the  field  strength  for  AC  and  point 
DC  dipoles  with  the  receiver  either  in  the  medium  or  above  the  medium.  In 
order  to  plot  the  special  case  of  a  finite  length  DC  dipole  with  the 
receiver  above  the  medium,  the  reduced  expression  had  to  be  derived.  These 
new  expressions,  in  addition  to  agreeing  with  Kraichman's  point  dipole 
expressions  which  are  valid  at  a  distance,  also  a*e  valid  with  the  receiver 
directly  above  the  source. 

ADM I N I STRAT I VE  I NFORMAT I ON 

This  project  was  supported  by  the  Annapolis  Laboratory  of  the  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center. 

INTRODUCTION 

Numerous  papers  have  been  published  to  give  various  approximate 
expressions  for  magnetic  field  strengths  in  various  configurations.  There 
is  a  growing  need  to  evaluate  the  magnetic  field  strengths  instantaneously 
and  interactively  so  that  curves  for  magnetic  fields  at  different  heights, 
different  distances,  different  frequencies,  and  different  orientations  in 
relationship  with  the  earth's  magnetic  field  can  he  produced  instantly. 

In  this  report,  the  author's  develope  a  set  of  computer  programs  to  evaluate 
the  magnetic  fields  generated  by  a  submerged  horizontal  electric  diode, 
dnly  the  quasi-stati  range  of  expressions  is  considered  throughout  the  report. 
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As  shown  in  Figure  1,  a  horizontal  electric  dipole  (HED)  of  finite 
length  L  angular  frequency  u,  is  located  at  depth  h  ( h  <  0 )  between  -L/2  and 

t 

L/2  parallel  to  the  x-axis  in  the  positive  x-direction.  The  plane  z= 0  separates 
the  upper  region  (z  >  0)  of  air  from  the  lower  region  of  the  conducting  medium 
with  conductivity  a  and  permeability  uQ.  The  observation  point  (x,y,z)  is 
either  in  the  air  or  in  the  medium.  For  simplicity,  Kraichman’s  aoproximate 

/  *  \  /  o  \ 

formulae^1'  and  Bannister  and  Dube's  modified  image  theory  results'  , 
instead  of  the  more  exact  numerical  results,  are  used.  A  collection  of  the 
reports  written  by  Bannister  and  his  co-worker  can  be  found  in  reference  3. 

The  expressions  for  a  finite  length  DC  dipole  are  derived  from  modified 
image  theory  results  as  a  limiting  case  of  an  AC  dipole.  Rectangular 
coordinates  are  used  throughout  the  report.  The  cylindrical  coordinate 
expressions  found  in  references  1  and  2,  are  charged  to  rectangular 
coordinates. 

The  computer  programs  are  written  in  77  Ansi  Fortran.  In  each  case  :  r 

an  AC  dipole,  the  computer  program  gives  the  moduli  of  each  of  the  three 

components,  the  moduli  of  the  component  in  the  direction  of  the  earth's 

magnetic  field,  the  real  part  of  this  projection,  and  the  phase  of  the 

projection.  In  all  of  the  AC  cases  the  frequency  is  assumed  to  be  1  Hz, 

the  dipole  current  is  50  amperes,  and  the  direction  or  the  earth's  magnetf 

field  is  in  the  negative  y  direction.  In  the  case  of  a  DC  dipole,  the 

computer  program  gives  the  three  components,  and  the  component  in  the 

direction  of  the  earth's  magnetic  field.  Curves  of  magnetic  fields  are 

(41 

produced  by  using  the  graphics  package'  '  developed  for  the  Textrcnix  4051 
at  the  U.  S.  Naval  Academy.  The  measuring  distances,  heights,  orientations, 
and  dipole  frequencies  can  be  specified  interactively  at  the  execution  time 
of  the  programs. 

Best  Available  Copy 
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The  programs  are  designed  so  as  to  generate  output  which  is  stored  in 
ta  files.  These  data  files  may  then  be  used  for  analysis,  or  they  may  Serve 
as  input  to  the  graphics  package  to  plot  the  curves.  This  approach  was  taken 
so  that  the  programs  can  be  used  even  if  plotting  capability  is  not  readily 
available.  A  complete  listing  of  all  the  programs  with  instructions  for  their 


use  can  be  found  in  Appendix  A. 


EQUATIONS 


Nine  different  configurations  of  field  strength  evaluations  will  be 
discussed.  The  first  six  cases  pertain  to  AC  sources.  Modified  image 
theory  results  are  used  in  the  first  three  of  these,  while  the  second  three 
describe  the  field  in  limited  ranges.  The  remaining  three  cases  refer  to  DC 
sources. 


(A)  Finite  length  AC  dipole  subsurface  to  air  propagation  (Modified 
Image  Theory). 

(2) 

Bannister  and  Dube  '  arrived  at  the  following  quasi-static  formulae 
by  employing  finitely  conducting  earth-image  theory  techniques. 

These  specific  expressions  are  in  simple  algebraic  form,  are  valid 
for  the  quasi-static  range,  and  can  be  found  on  page  in  of  reference  2. 


-  lle{ah  r  1  rd+z- 

(4>V  K12 


d+z-bh  z-bh. 


'22  (x^)S-y 


_1 _ .d+z-bh  z-bh. - 

fe)2+y2  K11  K21 


„  _  I  pyah  d+z-bh  x+2  x"2,  z-bh  f*T2  *"2,  AT2 

y  ‘  47f  y2+(d+z-bh)2  [K11  ‘  K12]  "  y2+(z-bh)2  Si  V  y2+(x+^ 

.d+z-bh  z-bh.. 

[  K11  "K22] 


X+y  X-y  ' 

f - - “1  + 

11/  1/  J 


Jj.erah  t  .1 _  r _ fL 

4,1  y2+(d+z-bh)2  K11 


x+k  x-t 


K12  y2+(z-bh)2  hi 


Bannister  and  Dube  chose  the  following  values  for  the  constants  a  and  b: 
a  =  0  and  b  =  1  for  R«6  and  (hj«6, 
a  =0.30  and  b  =  0.90  for  R/6  less  than  approximately  1, 
a  =  0.96  and  b  =  0.40  for  R/6  between  approximately  1  and  10, 
a  =  1.0  and  b  =  0  for  R>|3h|. 

The  name  of  the  program  used  to  predict  the  magnetic  field  strength  for 
this  case  is  ACF  (finite  length  AC  dipole  subsurface  to  air).  The  output  is 
shown  in  Figures  2,  3  and  4.  The  field  is  measured  along  the  path, 
x  =  -39.5m,  z  =  914.4m,  and  0<y^5000m.  Figure  2  depicts  the  absolute  values 
of  the  three  components  of  the  magnetic  field.  Figure  3  gives  the  modulus 
of  the  projection  of  the  magnetic  field  in  the  direction  of  the  earth's 
field  and  the  real  part  of  that  projection.  Figure  4  is  a  graph  of  the 


phase  of  the  projection. 
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yQy 


HJUi'f: 


(B)  Point  AC  dipole  subsurface  to  air  propagation  (Modified  Image 


Theory) , 


The  expressions  for  the  magnetic  field  subsurface  to  air 


propagation  due  to  a  point  AC  dipole  are  qiven,  in  cyclindrical  coordinates 
on  page  6  of  reference  2.  As  pointed  out  in  reference  2,  these  expressions 
can  be  derived  from  the  corresponding  finite  length  DC  expressions  (1),  (2), 
and  (3),  as  the  measurement  distance  becomes  much  greater  than  the  source 
length  L.  In  the  interest  of  consistency,  the  expressions  given  in 
reference  2  have  been  transformed  into  rectangular  coordinates. 

,,  „  Px.v  _  yah  fd+z-bh  z-bh  A  2  ,d+z-bh  z-bhv, 

"x  ?  e  l  -3  “  "3  +  v  V  ~  ~V  /  i 


x  4,p2 


k3  p2  K2 


u  _  P  „  yah  r /d+z-bh  z-bh^  y2  ^  y2-x2  ;d+z-bh  z-bh^ 
'■y  zhre  7  T~>  75  +  ~7a —  \ — i? - i~~~S  I 


i/3  „3  '  p2  p4  '  K,  K,  n 

h  h  2i  (5) 


>z  -  4V7 


%  vah  (_1_  _  u 
4tt  *3  y3 

h  K2 


ihe  constants  a  and  b  are  again  loosely  defined  as  in  Case  A.  The  name  of 
the  program  which  produces  the  output  for  this  case  is  ACP  (point  AC  dipole 
subsurface  to  air).  Sample  outputs  of  ACP  are  shown  in  Figures  5,  6  and  7. 

All  the  constants,  as  well  as  the  three  output  examples,  are  the  same  as  the 
previous  representation  for  a  finite  dipole. 


(C)  Point  AC  dipole  subsurface  to  subsurface  propagation  (Modified 
Image  Theory). 

The  expressions  for  the  magnetic  field  subsurface  to  subsurface 
propagation  due  to  a  point  AC  dipole  are  given,  in  cylindrical  coordinates, 
on  pages  13  and  14  of  reference  2.  In  rectangular  coordinates,  they  become: 


H.  -  {iih  (i+YR2)e-^R2  ♦  t  2b(z^L  + 


4  Tip  RJ 


k4p‘ 


o  K- 


(7) 


H  =  JL  {.  £jl 

Hy  «,  '  R  3 


(lnR,)e-YRl  -  (l*YR,)e-YR2  *  e  a(z+h)  [^^1  (y2-x2) 

1  p^R^  l/  - 


V 


+  b(z+h)(y2-x2)  d-b(z+n)  A } 

4  ,3  2J/ 

k4  p 


K3P 


(8) 


=  {^zr  {1+yRi} 


-  R 


-  <L_2  (1+Y»  )  +  eva(^),  1  .  1 


(9) 


The  constants  a  and  b  are  determined  by: 
a  =  0  and  b  =  1  for  R^/c^-U, 

a  =  0.4  and  b  =  0.95  for  R^/6  less  than  approximately  1, 
a  =  0.96  and  b  =  0.4  for  between  approximately  1  and  10, 
a  =  1  and  b  =  0  for  p>3|z  +  hj. 

The  name  of  the  program  is  ACPS  (point  AC  dipole  subsurface  to  subsurface). 
The  sample  outputs  for  this  case  are  shown  in  Figures  8,  9,  and  10.  Since  this 
example  is  for  the  case  where  both  the  dipole  and  the  sensor  are  subsurface,  the 
value  of  z  is  now  z  =  -5m.  The  path  y  is  once  again  chosen  as  0<y; 5000m. 
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POLE  SUBSUPFAC 


Yj  HI— SOLID? HR-DOT 


(D)  Point  AC  dipole  subsurface  to  air  propagation  in  the  range 
|yRJ»1  and  R>>  jhj . 

The  next  three  sets  of  equations  (D),  (E),  and  (F)  are  valid  only 
in  restricted  ranges. 

The  expressions  for  the  magnetic  field  subsurface  to  air  propagation 
in  the  range  { yR | »1  and  R»|h|  are  given,  in  cylindrical  coordinates  on 
page  3-22  of  reference  1.  In  rectangular  coordinates,  they  become: 


H  =  ip**.  s^L  { iji) 

*  ^  7? 1  r2’ 


P  eyh  /o  .2  2  3z2y2 


"y  =  ^  R3p 


/  O  C.  2 

3.2  ^2y  "X  _  r2 


=  (1+Y2.5j!> 

2tiy  R°  R^ 


(10) 


(11) 


(12) 


The  name  of  the  program  for  case  (D)  in  ACPI  (Point  AC  Dipole 
Subsurface  to  Air).  Figures  11,  12,  and  13  are  representative  outputs  of 
this  program. 
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I  NT  AC  DIFOLE  SPECIAL  PhHGE  #1 


{£)  Point  AC  dipole  subsurface  to  air  propagation  in  the  range 


J 

I 

|yp|»1.  P>>  i  h  I ,  P»  |  z  |  • 

I  The  expressions  for  the  magnetic  field  subsurface  to  air  propagation 

in  the  range  | yp  1  >>  1 ,  p»|h|  and  p»|z|  are  given,  in  cylindrical  coordinates, 
on  page  3-24  of  reference  1.  The  expressions  are  valid  only  in  the  vicinity 
of  the  z  axis.  In  rectangular  coordinates,  they  become: 


3Pxy  eyh 

2iryp 

(13) 

PpY^  2  ? 

*2-5  (2y -x2) 

27TYP 

(14) 

eYW) 

2ir/pb 

(15) 

The  program  for  these  equations  is  named  ACP2  (Point  AC  Dipole 
Subsurface  to  Air).  Figures  14,  15,  and  16  are  included  as  examples.  The 
graphs  are  plotted  along  the  path  x  -  152.4m,  z  =  152.4m,  and  0<y55000m. 


21 

I 


POINT  AC  DIPOLE  SPECIAL  RANGE  #2 


,  HX-SOL 1 D , HY-DhSH? HZ-DOT 


DIPOLE  SPECIAL  RANGE  #2 


SOLID, HR-DOT 


POINT  hC  DIPOLE  SPECIhL  RANGE  #2 


(F)  Point  AC  dipole  subsurface  to  subsurface  propagation  in  the  range 
!yp!»1,  p»'|h|,  and  p»|z|. 

The  range  of  this  case  is  identical  to  that  of  case  (E).  The 
difference  lies  in  the  fact  that  the  observation  point  is  in  the  lower  half 
medium  as  opposed  to  the  upper  half* 

The  expressions  for  the  magnetic  field  subsurface  to  subsurface 
propagation  in  the  range  |yp|»1,  p»|hj,  and  p»jzj  are  given,  in 
cylindrical  coordinates,  on  page  3-24  of  reference  1.  In  rectangular 
coordinates,  they  become: 

H  -  3Pxyp^z**^ 


V 


0  5 

2~yp 

(16) 

pAY(z+h)  2  2 

5-  (2/-X2) 

2ttyp 

(17) 

3Py  eY(h+z) 

2  5  e 

2ity  p 

(IS) 

set  of  equations  has 

been  named  ACP2S  (Point 

AC  Dipole  Subsurface  to  Subsurface).  The  sample  outputs  are  shown  in 
Figures  17,  18,  and  19.  The  path  of  measurement  is  x  =  152.4m,  z  =  -5m, 
and  0<y^500Cm. 
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POINT  AC  DIPOLE  SPECIAL  RANGE  #2  SUBSURFACE 


(G)  Point  DC  dipole  subsurface  to  air  propagation. 


The  remaining  three  programs  apply  to  DC  sources  only.  These 
programs  have  been  used  extensively  by  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center,  Annapolis,  Laboratory. 

The  expressions  for  the  magnetic  field  subsurface  to  air  propagation 
due  to  a  point  DC  dipole  are  given,  in  cylindrical  coordinates,  on  page  3-4 
of  reference  1.  In  rectangular  coordinates,  they  become: 


H*  -  f-i  fir +  h  <sf  - 

4 -no  Rj  o  1 


V 


_P_  r(z-h)y"  x2-y2  ,2-h 
*  1  _  3  '  2  lR. 


4-0 


- 1)] 


{19} 

(20) 


H 


z 


(21) 


The  program  for  this  set  of  equations  is  called  DC?  (Point  DC  Dipole 
Subsurface  to  Air  Propagation).  The  sample  graphs  for  this  case  are  labeled 
Figures  20  and  21.  The  path  for  Figure  20  is  x  =  304.8m,  z  =  914.4m,  and 
-2500SyS2500.  Since  the  equation  has  a  singularity  at  the  origin,  x  was  chose 
to  be  0.01m.  The  variables  z  and  y  are  the  same  for  both  figures. 
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(H)  Finite  length  DC  dipole  subsurface  to  air  propagation. 
The  magnetic  field  expressions  ( 1) ,  (2),  and  (3)  for  a 
finite  length  AC  dipole  can  be  used  to  derive  the  corresponding  DC 
expressions.  As  d^,  e'a^+l,  and  b  =  1,  we  have: 


d+z-bh 

^12 


d+z-bh  , 
kll 


-0, 


K12 


->0. 


and  equations  (1),  (2),  and  (3)  reduce  to: 

(22) 

(23) 

(2^) 

Equations  (22),  (23),  and  (24)  give  the  magnetic  field  components 
for  a  finite  length  DC  dipole.  When  the  finite  length  L  is  much  less  than 
the  measurement  distance,  we  have: 


H  = 


H  = 


ll  r _ I _ 

4u  L,  L , 2 ,  2 


<i-kt)  - 


1 


{x-^+y"  ^22  (x+|)2+y2 


(1  -  £*)] 

72 


J_  r  z-h 
4ti  ‘ .  2" 


x+^-  x~Y 

7  t '  '  ;.2" 


y  +(z-hr  n21  n22  y4+(x+|) 


X+?  7  h 

,-(i-r)! 
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H_  =  - 


iy 


4t?  y2+(z-h)d  n21 


x^  x— Sj- 

2  (icr " 


(1-  —o')* 

2R, 


1 

K22 


1 

1 


(1- 


_xL  . 

?  >. » 

2R^ 

a 


4  (i  +  4). 


(x-2)‘ 


I 


and  IL  -*  P. 


Consequently,  equations  (22),  (23),  and  (24)  reduce,  respectively,  to  the 
point  dipole  equations  (19),  (20),  and  (21). 

The  name  of  the  program  for  this  set  of  equations  is  DCF  (Finite 
Length  DC  Dipole  Subsurface  to  Air). 

Figure  22  shows  the  magnetic  fields  due  to  a  DC  dipole  of  length 
50m,  located  at  a  depth  of  76.2m,  for  the  path  z  =  914.4m,  x  =  304.8m,  and 
-2500m5y-2500m.  The  earth's  magnetic  field  direction  in  the  xy-plane  forms 
an  angle  of  1.5ir  radium  relative  to  the  x  axis.  Figure  23  gives  the 

corresponding  curve  for  a  point  DC  dipole.  Figure  22  and  Figure  23  show  that 
along  this  path,  the  finite  length  dipole  results  are  in  very  good  agreement 
with  corresponding  point  dipole  results.  The  point  dipole  equations  (19),  (20), 
and  (21)  have  a  singularity  at  the  origin.  Therefore,  they  are  not  valid  near 
the  z  axis.  Figures  (21)  and  (23)  point  out  the  discrepancy  encountered  when 
the  observation  path  is  near  the  z  axis.  Both  figures  are  measured  at  the 
same  height,  z  =  914.4m,  and  y  is  bounded  by,  -2500<y<2500.  Even  though  the 
values  for  x  are  very  close  to  each  other,  0.00m  and  0.01m,  the  magnetic 
field  strengths  appear  to  be  vastly  different  near  the  z  axis. 

The  finite  dipole  equations  (22),  (23),  and  (24)  predict  magnetic 
fields  close  to  that  predicted  by  the  point  dipole  equations  when  the 
observation  points  are  away  from  the  origin.  In  addition,  they  do  not  have 
a  singularity  at  the  origin. 
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FINITE  LENGTH  DC  DIPOLE,  BETA*1 


580  -500  500  1500  2500 

■SOLID?  HY-DtiSH?  H2-CHAIN*  HT-DOTTED?  2-914. 4, Xs8. 0 


(I)  Point  DC  dipole  subsurface  to  subsurface  propagation. 

The  expressions  for  the  magnetic  field  subsurface  to  subsurface 
propagation  due  to  a  point  DC  dipole  are  given,  in  cylindrical  coordinates, 
on  page  3-4  of  reference  1.  In  rectangular  coordinates,  t.iey  become: 


H  =  +  2_  ,z+h  , 

x  .  2  lD  3  2  [R0  i}‘ 

p  R~  p  2 


(25) 


,,  -  P  ,z-h  ,  x2(z+h)  x2-y2  fz+h  n , 
Hy  4¥(7T  TT— +  4"  (RT  +  1}j 
J  -  p  R2  0  z 


(26) 


H 


z 


Py 

4ttR13 


(27 


The  name  of  the  program  is  DCPS  (Point  DC  Dipole  Subsurface  to 
Subsurface).  The  appropriate  graph  is  given  in  Figure  24.  The  path  of 
measurement  is  z  =  -5m,  x  =  -39.5m,  and  0  S  y  S  SOOOn. 
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DISCUSSION 


The  computer  programs  can  be  used  to  find  the  magnetic  fields  along 
certain  chosen  paths  in  the  upper  half  or  in  the  lower  half  medium.  The 
corresponding  magnetic  curves  for  each  path  can  then  be  plotted.  The  sys 
which  predict  the  results  for  DC  sources  have  been  utilized  extensively 
by  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center,  Annapo 
Maryland. 


APPENDIX  I 


Operating  Instructions  and  Program  Listings 
OPERATING  INSTRUCTIONS 

The  programs  are  intended  to  be  run  in  an  interactive  mode.  Data  which 
is  required  for  each  individual  program  is  requested  by  the  program  at  run 
time.  Depending  on  the  particular  program  which  is  being  run,  the  user  will 
be  required  to  furnish  the  values  of  some  combination  of  the  following 
parameters: 

(a)  The  height  of  the  measurement  path,  that  is,  the  value  of  2. 

(b)  The  values  for  the  constants  a  and  b,  required  when  using  Modified 
Image  Theory. 

(c)  The  range  of  the  y  value  of  the  measurement  path  and  the  increment 
along  that  path. 

(d)  The  angle  of  the  earth's  magnetic  field  direction.  The  angle  is 
assumed  to  be  in  K  -  radians,  but  the  user  need  only  provide  the 
value  for  K. 

The  output  of  the  programs  is  stored  in  data  files.  Each  AC  program 
writes  into  three  files;  one  for  Hx,  and  Hz;  one  for  Ht  and  Hr;  and  the 
third  is  used  to  store  the  phase.  Each  DC  program  requires  only  one  output 
file.  H  ,  H  ,  H  and  H.  are  all  stored  in  this  file.  These  files  may  now 

X  jr  Z  U 

provide  input  to  TEKGRAF2,  a  graphics  plotting  package  developed  at  the 
United  States  Naval  Academy.  Reference  3  provides  operating  instructions 
for  this  package. 
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ACF 

♦  MODIFIED  IMAGE  THEORY  FOR  FINITE  LENGTH  AC  SUB  TO  AIR 

*  NUSC  REPORT  5647 

COMPLEX  FUNCTION  HX<FREQ»X»Y,Z) 

COMPLEX  ZDBrZDBY »K11*K12 
COMPLEX  COEFF  r GAMMA 
REAL  K21 f K22 
REAL  LENGTH 

COMMON  DEPTH  *  LENGTH  »  CURENT  »  ACON , BCON  »  PI 
GAMMA  =  CSQRT<CMPLX<0.0fFREQ*<PI**2)*32.E-7> ) 

XLP  =  X  +  LENGTH/2 

XLM  =  X  -  LENGTH/2 

ZB  =  Z  -BCON*DEPTH 

ZDB  =  Z  +  2/GAMMA  -  BCON*DEPTH 

XLPY  =  XLP**2  +  Y**2 

XLMY  =  XLM**2  +  Y**2 

ZDBY  =  ZDB#*2  +  Y**2 

ZBY  =  ZB**2  +  Y**2 

Kll  =  CSQRT <XLP**2  +  ZDBY) 

K12  =  CSQRT <XLM**2  +  ZDBY) 

K21  =  SQRT<XLP**2  +  ZBY) 

K22  =  SQRT(XLM**2  +  ZBY) 

COEFF  =  Y*CURENT  *  CEXP(GAMMA*ACON*DEPTH> 

1  HX  =  100.*C0EFF*(<ZDB/K12  -  ZB/K22)/XLMY  -  (ZDB/K11  - 
RETURN 

END 

COMPLEX  FUNCTION  HY<FREQ»X*Y»Z) 

COMPLEX  ZDB » ZDBY »K11»K12 
COMPLEX  HY1 *HY3»HY4»GAMMA» YCOEFF 
REAL  K21»K22 
REAL  LENGTH 

COMMON  DEPTH  »  LENGTH  »  CURENT  *  ACON  »  BCON  r  F  I 
GAMMA  =  CSGRT<CMPLX<O,0fFREQ*<PI**2)*32.E-7) ) 

XLP  =  X  +  LENGTH/2 

XLM  =  X  -  LENGTH/2 

ZB  =  Z  -BCQN*DEPTH 

ZDB  =  Z  +  2/GAMMA  -  BCON*DEPTH 

XLPY  =  XLP**2  +  Y**2 

XLMY  =  XLM**2  +  Y**2 

ZDBY  =  ZDB**2  +  Y**2 

ZBY  =  ZB**2  +  Y**2 

Kll  =  CSQRT  <XLP$*2  +  ZDBY) 

K12  =  CSQRT <XLM**2  +  ZDBY) 

K21  =  SORT (XLP*#2  +  ZBY) 

K22  =  SQRT<XLM**2  +  ZBY) 

YCOEFF  =  CURENT  *  CEXP< GAMMA*ACON*DEPTH ) 

2  HY1  =  (XLP/K11  -  XLM/K12)*ZDB/ZDBY 
HY2  =  (XLP/K21  -  XLM/K22)*ZB/ZBY 
HY3  =  (ZDB/K11  -  ZB/K21 XXLP/XLPY 
HY4  =  (ZDB/K12  -  ZB/K22)*XLM/XLMY 
HY  =  100 ♦ #YCOEFF* ( HY1-HY2+HY3-HY4 ) 

RETURN 

END 
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COMPLEX  FUNCTION  HZ (FREQf X »  Y  *  Z ) 

COMPLEX  ZDBfZDBY»KllrK12 
COMPLEX  GAMMA  *  COEFF 
REAL  K21 t K22 
REAL  LENGTH 

COMMON  DEPTH * LENGTH  *  C'JRENT  *  ACON  »  BCON » PI 
GAMMA  =  CSQRT ( CMPLX  <0.0*  FREQ* ( PI#*2 ) *32 . E-7 ) ) 

XLP  =  X  +  LENGTH/2 

XLM  =  X  -  LENGTH/2 

ZB  =  Z  -BCON*D£PTH 

ZDB  =  Z  +  2/GAMMA  -  BCON*DEPTH 

ZDBY  =  ZDB**2  +  Y**2 

ZBY  =  ZB**2  +  Y**2 

Kll  =  CSGRT(XLP**2  r  ZDBY) 

K12  =  CSGRT(XLM**2  +  ZDBY) 

K21  =  SGRT(XLP**2  +  ZBY) 

K22  =  SQRT(XLM**2  +  ZBY) 

COEFF  =  Y*CUR£NT  *  CEXP(GAMMA*ACON*DEPTH) 

3  HZ  =  100.*C0EFF*< (XLP/K11  -XLM/K12 ) /ZDBY  -  ( XLP/K21  -  XLM/K22)/Zc. 

4  RE  TURN 
END 

COMPLEX  FUNCTION  HT (FREQ* X» Y »ZrBETA) 

COMPLEX  HXrHYfHZ 
SBETA  =  SIN (BETA) 

CBETA  =  COS (BETA) 

HT  =  HX(FREQ.XrYfZ)*CBETA/l .76  ~0.823*HZ(FREGf Xf Y,Z) 

HT  =  HT  +  HY (FREQf X*Y  »Z)#SBETA/1 .76 

RETURN 

END 

FUNCTION  HP(FREQrXfYfZfBETA) 

COMPLEX  HXrHYfHZrHT 

HP  =  ATAN( AIHAG(HT (FREQfXf Y  fZrBETA) ) /REAL  <  HT ( FREQ  rX»Y»Z»BEfA>) > 

RETURN 

END 

FUNCTION  HR(FREQrXf Y f Zf BETA) 

COMPLEX  HXfHY,HZ 
SBETA  =  SIN( BETA ) 

CBETA  =  CCS (BETA) 

HR  =  R£AL(HX(FREQ»XfYfZ) >*CBETA/1.76  -  0 . 823*REAL ( HZ ( FREQ r  X  f  Y  f  Z  >  ) 

HR  =  HR  +  REAL ( HY ( FREQ r X *  Y *7  >  > *SBETA/1 .76 

RETURN 

END 

COMMON  DEPTH  *  LENGTH  f  CURENT » ACON  r  BCON  f  PI 

COMPLEX  HXfHYyHZf HT 

REAL  LENGTH 

DEPTH  =  -76.2 

LENGTH  =  50. 

CURENT  =  1.00 
PI  =  3.1415? 

X  =  -39,5 

OPEN  (If' ACFFD ' f  ACCE3S=  *  ASCII  *  ) 

OPEN  ( 2  f *  ACFPFD ' r  ACCESS=  *  ASC II*) 
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OPEN  <  3  *  *  ACFTFD  *  * ACCESS=  *  ASC II") 

FREQ  =  1 

PRINT*  ’WHAT  IS  THE  HEIGHT  IN  METERS’ 

INPUT*  Z 

PRINT *  ’WHAT  IS  THE  LIMIT  FOR  Y  AND  ITS  INCREMENT’ 
INPUT  *  LIMIT  »  INC 

PRINT *  ’WHAT  ARE  THE  VALUES  FOR  A?  AND  FOR  B’ 

INPUT *  ACONfBCON 

PRINT *  ’WHAT  IS  THE  ANGLE  FOR  EARTH  MAGNETIC  FIELD’ 

INPUT *  ANGLE 

BETA  =  ANGLE*PI 

DO  10  Y  =  0*01 * LIMIT * INC 

AHX  =  CABS ( HX ( FREQ  *X*Y*Z)) 

10  WRITE  (1*100)  Y*AHX 
WRITE  ( 1 *200) 

DO  20  Y  =  0 « 01 r LIMIT* INC 
AHY  =  CABS  <  HY  <  FREQ  *X*Y*Z)> 

20  WRITE  < if  100)  Yt AHY 
WRITE  <1*200) 

DO  30  Y  =  0.01»LIMIT*INC 
AHZ  =  CABS ( HZ ( FREQ  *X*Y*Z>) 

30  WRITE  (1*100)  Y*  AHZ 
WRITE  (1*200) 

DO  40  Y  =  0*01  * LIMIT  *  INC 

AHT  =  CABS(HT (FREQ*X*Y*Z*BETA) ) 

40  WRITE  (3*100)  Y*  AHT 
WRITE  (3*200) 

DO  50  Y  =  0 .01  * LIMIT *  INC 
50  WRITE  (3*100)  Y*  HR(FREQ*X* Y*Z*BETA) 

WRITE  (3*200) 

DO  60  Y  =  0*01* LIMIT* INC 
60  WRITE  (2*100)  Y*  HP(FREQ*X*Y*Z*BETA) 

WRITE  (2*200) 

CLOSE  (1) 

CLOSE  (2) 

CLOSE  (3) 

100  FORMAT  <1PF10.3*1H*  *1PE10.3) 

200  FORMAT  <11H1.E37*1.E37) 

END 
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MODIFIED  IMAGE  THEORY  FOR  INFINITESIMAL  AC  HF.D 
SUBSURFACE  TO  AIR  PROPAGATION 
BANNISTER  REPORT  5647 
COMPLEX  FUNCTION  HX(FREQrXrY rZ> 

COMPLEX  ZDBrK2rHXlrHX2rGAMMA 
COMPLEX  COEFF 
REAL  Kl-PI 

COMMON  DEPTH r ACON t BOON r P f PI 

GAMMA  -  CSQRT (CMPLX <0. Or  FREQ* <PI**2 >*32. E~7> ) 

ZB  -•=  7.  -BCON*DEPTH 

ZDB  =  Z  +  2/GAMMA  -BCON&DEPTH 

COEFF  =  P  *  CEXP(GAMMA*ACUN*DFprH)/(4*PI> 

RHO  =  SGRT(X**2  4-  Y**2> 

K 1  a  SORT ( RHO**2  +  ZB**2> 

K2  =  CSQRT <RHG**2  +  ZDB**2> 

HX1  =  ZDB/(K2**3>  -  ZB/(K1**3> 

HX2  =  (ZDB/K2  -  ZB/K1 ) *2/ ( RH0**2 ) 

HX  =  400.*PT*C0EFF*X*Y*<HXt  +  HX2 ) / ( RH0**2 ) 
RETURN 


END 

COMPLEX 

COMPLEX 

COMPLEX 

COMPLEX 


FUNCTION  HY(FRFOfXfYrZ) 
ZDB rK2.HYlfHY2f GAMMA 


HY1 » HY2  r HX1 f  HX2  r  GAMMA 
COEFF 
REAL  Kir  PI 

COMMON  DEPTH r ACON r BCON » P r PI 

GAMMA  =  CSQRT ( CMPLX < 0. 0 r FREQ* <PJ**2> *32. E -7 > ) 

ZB  =  Z  -BCON*DEPTH 

ZDB  =  Z  F  2/GAMMA  --  BCONftDEPTH 

RHO  =  SGRT(X**2  +  Y**2) 

K 1  =  SORT ( RH0**2  +  ZB**2> 

K2  =  CSQRT  ( RHQ**2  +  ZDB**2> 

COEFF  =  P  4.  CEXP<GAMMA*AC0N*DEPTH)/<4*PI> 

HYi  =  < ZDB/ ( K2*K2*K2 )  ~  7B/(K1*K1*K1 ) >*X*Y/(RH0**2> 
HY2  =  (ZDB/K2  -  ZB/Ki)«<YS*2  -  X#*2>/<PHn**4> 

HY  =■  400.*PI#C0EFF*<HY1+HY2> 


END 

COMPLEX  FUNCTION  HZ (FREQ *X  r  Yr  Z  > 
COMPLEX  ZDB rK2* GAMMA 


REAL  KlrPl 

COMMON  DEP  TH  r  ACON  r  BCON  r  P  ?  P I 

GAMMA  =  CSQRT ( CMPLX(0. Or FREQ* <PI**2> *32. E~7>> 

ZB  a  z  -BCONXcBEPTH 

ZDB  =  Z  +  2/GAMMA  -  BCON*DEPTH 

RHO  a  SORT ( X**2  +  Y**2> 

K1  =  SQRT ( RM0**2  +  ZB**2> 

\2  ~  CSPRT ( RH()**2  +  ZDB**2> 

COEFF  a  P  *  CEXP(GAMMA*ACON*DEFTH)/<4*PI ) 

HZ  a  400.*PI*Y*C0EFF*CI/<K1**3>  -  2/<K2**3>> 

RETURN 

END 

COMPLEX  FUNCTION  HT ( FREQ r X r Y r Z r BETA > 
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COMPLEX  HX*HY*HZ 
SBETA  =  SIN( BETA) 

CBETA  =  COS (BETA) 

HT  =  HX ( FREQ *  X  *  Y  *  Z ) #CBETA/1 ♦ 76  i  HY (FREQ > X  * Y » Z > *SBETA/1 . 76 
HT=  HT-0 . 823#HZ  <  FREQ*  X *  Y  *  Z  ) 

RETURN 

END 

FUNCTION  HP ( FREQ *  X  *  Y  *  Z  *  BETA ) 

COMPLEX  HX*HY*HZ*HT 

HP  =  A" AN< AIMAG(HT (FREQ*X  * Y  *  Z  * BETA) ) /REAL(HT ( FREQ*  X •  Y  *  Z  *BETA  > ) ) 

RETURN 

END 

FUNCTION  HR  <  FREQ *  X  *  Y  *  Z  * BET  A ) 

COMPLEX  HX*HY*HZ 
SBETA  =  SIN(BETA) 

CBETA  =  COS (BETA) 

HR  =  REAL < HX( FREQ >X»Y  *Z) )*CBET A/ 1 .76 
HR  =  HR+REAL(HY(FREQ*X»Y*Z) )#SBETA/1 .76 
HR  =  HR  --  0 * 823#REAL ( HZ ( FREQ *  X  *  Y  *  Z )  ) 

RETURN 

END 

COMMON  DEPTH *ACQN»BCON»P* PI 
COMPLEX  HXrHY  *HZ»HT 
X=-39 . 5 
DEPTH  =  -76 « 2 
PI  =  3.14159 
P  =  50.0 

OPEN  ( 1  * 'ACPFD' *  ACCESS” 'ASCI I = ) 

OPEN  ( 2 * * ACPPFD  * » ACCESS* ' ASCI I  *  > 

OPEN  ( 3 * ■ ACPTFD  *  r  ACCESS* 'ASCII  *  ) 

FREQ=1 ♦ 0 

PRINT *  'WHAT  IS  THE  HEIGHT  IN  METERS' 

INPUT  *  Z 

LIMIT=5000 

INC=100 

PRINT *  'WHAT  ARE  THE  VALUES  FOR  A?  AND  FOR  B* 

INPUT *  ACGN?BCGN 

PRINT *  'WHAT  IS  THE  ANGLE  OF  THE  EARTH  MAGNETIC  FIELD' 

INPUT  *  ANGLE 

BETA  =  ANGLE*PI 

DO  10  Y  =  0.01 r LIMIT » INC 

AHX  =  CABS ( HX ( FREQ  *  X  *  Y  *  Z )  > 

10  WRITE  (1*100)  Y*  AHX 
WRITE  ( 1  * 200 ) 

DO  20  Y  =  0«01*LIMIT* INC 
AHY  =  CABS(HY(FREQ*X*Y*Z) ) 

20  WRITE  (1*100)  Y »  AHY 
WRITE  ( 1*200) 

DO  30  Y  =  0.01*LIMIT* INC 
AHZ  =  CABS(HZ(FREQ*X* Y  *Z) ) 

30  WRITE  (1*100)  Y*  AHZ 
WRITE  (1*200) 
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*  MODIFIED  IMAGE  THEORY  FOR  POINT  AC  SUB  TO  SUB 

*  NUSC  REPORT  564/ 

COMPLEX  FUNCTION  HX ( FREQr X  »  Y  »  Z  > 

COMPLEX  HX1»MX2 

COMPLEX  EXP2  r EXPZ , GAMMA  »  DZB , KFOUR 
REAL  K THREE 


COMMON  DEPTH  r  A CON »  BOON r P  r P I 

GAMMA  ~  CSQRT<CMPLX(O.OrFREQ*<PI** 2  *32,E~7)> 

ZHP  =  Z  4  DEPTH 

DZB  »  2/GAMMA  -  DC ON* ZHP 

RHO  =  SORT ( X**2  +  Y**2> 

RTUO  =  SORT  ( RHQ#*2  +  ZHP**2> 

EXP2  =  CEXP ( -GAMMA*RTWQ > 

EXPZ  “  CEXP<GA«MA*ACON*ZHP> 

K THREE  =  SGRT<ftM0»*2  +  ( BCONftZHP > **2 > 

Kr OUR  =  CSGRT<RHO**2  +  BZB**2) 

HX 1  =  ZHP* (i  r  GArtMA*RTUO>*EXP2/(RTWO**3> 

H>  2  -  2*BZB/ ( KFOUR* ( RHO**2 ) >  +  2*BC0N*ZHP/ (KTHREE* (RH0**2 ) > 
MX  2  =  HX2  +  BZ8/(KF0UR**3> 

MX  =  100.*P#X*Y»(HX1  *  EXFZ*HX2>/(RH0**2> 

RETURN 


END 

COMPLEX  FUNCT I ON  HY ( FREQ ,XrYrZ) 

COMPLEX  HY  1  -t  HY2  *  HY3 » EXP1 
COMPLEX  EXP2 / EXPZ 7  GAMMA  r DZB  r KFOUR 
REAL  KTHREE 

COMMON  DEPTH ; ACON t BOON r P r PI 

GAMMA  s  CSGRT (CKFLX(0. Or  FREQ* (PI**2> *32 . E-7 ) ) 

ZHM  =  Z  -  DEPTH 

ZHP  =  7  +  DEPTH 

DZB  =  2/GAMMA  -  DC ON* 7 HP 

RHO  --  SQRT<X#*2  ■}•  Y**2> 

RONE  «  SORT  <  RHO* *2  +  ZHM**2> 

RTWO  “  SORT <RHO**2  +  ZHP**2> 

EXP1  ==  CEXP(~GAMMA*RQNE> 

EXP2  -  CEXP<-GAMMA*RTWO> 

EXPZ  =  CEXP  <  GAMMA# ACON*ZHP ) 

KTHREE  =  SQRf<?::-:0**2  +  <  BCON*ZHP )  **2  > 

KFOUR  =  CS9RY<SH0**2  +  DZ8**2> 

HY 1  =  -ZHM* < i  >  GAHNA*R0NE)#EXP1/(R0NE**3> 

HY2  =  “ZHP*(X**2)#<1  +  GAMMA*RTWO)*EXP2/( ( RH0**2  >  *  ( RT«JO**3  > ) 
HY3  =  DZB#(Y**2  -  X**2)/(KF0UR*<RH0**4> > 

HY3  =  HY3  t  8C0N*ZHP*(Y**2  -  X**2> /(KTHREE* <RH 0**4 >  ) 

HY3  =  £XPZ*(HY3  T  ( Y**2 > *DZB/ ( ( KF0UR**3 ) * < RHO* Y2> ) ) 

HY  =  1 00 ♦ ( HY 1  +  HY2  +  HY3 ) 

RETURN 

END 

COMPLEX  FUNCTION  HZ ( FREQ y X c Y r Z ) 


COMPLEX  HZ  i  *  B2'2  v  HZ3 1 EXP1 
COMPLEX  EXP2 y EXPZ r GAMMA r  DZB  r  KFOUR 
REAL  KTHREE 

COMMON  DEPTH r ACON r BCON r  P  y  PI 
07/21/80 


G°^ 

14  ♦ 06  J 40 


A7 


ACPS 


(Pa<3e  2) 


GAMMA  =  C3QRT(CMPLX(0.0?FREQ*(PI**2>*32.E-7> ) 

ZHM  =  Z  -  DEPTH 

ZHP  =  Z  +  DEPTH 

DZB  =  2/GAMMA  -  BCON*ZHP 

RHO  =  SORT  ( X**2  4  Y**"" 

RONE  =  SORT  ( RH0**2  -I  z.HM**2) 

RTWO  =  SORT ( RH0**2  +  ZHP**2) 

EXPi  =  CEXP ( -GAMMA*RONE ) 

EXP2  =  CEXP < ~GAMMA*RTUO ) 

EXPZ  =  CEXP <  GAMMA*ACQN*ZHP ) 

MHREE  =  SORT ( RH0**2  +  (BC0N*ZHP>**2 > 

KFUUP  =  CSQRT ( RHO**2  4-  DZB**2) 

HZ1  -  EXPl* ( 1  +  GAMMA*RONE ) / ( R0NE**3 ) 

HZ2  =  EXP2*  < 1  +  GAMMA*RTWO  )  /  ( RTW0**3  > 

HZ3  =  EXPZ*  ( t/(,KTHREE**3  )  -  i /  < KF0UR**3 )  ) 

HZ  =  lOO.*P*Y*<HZl  -  HZ2  +  HZ3)/RH0 
4  RETURN 
END 

COMPLEX  FUNCTION  HT ( FREQ ? X ? Y ? Z s- BETA ) 

COMPLEX  HXeHY.HZ 
SUET A  =  SIN (BETA) 

CBETA  =  COS (BETA) 

HT  =  HX(FREQ?X?Y  »Z)*CBETA/1 .76  -  0 . 823*HZ ( FREQ t X , Y , Z ) 

HT  a  HT  +  HY(FREQ?X?Y ?Z)*SBETA/1 .76 

RETURN 

END 

FUNCTION  HP(FREQ?X? Y  rZ?BETA) 

COMPLEX  HX  ?  MY  ?  HZ  ? HT 

HP  =  AT  AN ( A1MAG ( HT ( FREQ  ?  X  ?  Y  r Z  r  BET  A ) ) /REAL ( HT ( FREQ » X  >  Y » Z  f  BETA 

RETURN 

END 

FUNCT  ION  HR ( F  REQ  >  X  ?  Y*  7  r BETA  > 

COMPLEX  H\ • HY  » HZ 
SBETA  -  SIN  <  BETA ) 

CBETA  =  COS (BETA) 

HR  =  REAL(HX(FREQrXrYrZ) )*CBETA/1 .76  -  0 . 823*P.AL ( HZ ( FREQ ? X ? 

HR  =  HR  +  REAL ( HY ( FREQ  ?  X  ?  Y  ?  Z ) ) *SBET A/ 1 .76 

RETURN 

END 

COMMON  DEPTH t ACON  r BOON  ? P  r  P I 
COMPLEX  HXrHY  ?HZ»HT 
X=76 ♦ 2 

DEPTH  =  -76.2 
PI  =  3.14159 
P  =  50.0 

OPEN  < 1 r * ACPSFD ■ ? ACCESS= *  ASC II*) 

OPEN  .2* “mCPSPFD* ?ACCES3=*ASCII * ) 

OPEN ( 3 ? -ACPSTFD* »ACCESS=’ ASCII'  ) 

FREQ=1 

PRINT r  5 WHAT  IS  THE  HEIGHT  IN  METERS* 

INPUT  r  Z 

PRINT?  "WHAT  IS  THE  UPPER  LIMIT  FOR  Y  AND  ITS  INCREMENT* 
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INPUT ,  LIMIT, INC 

PRINT ,  “WHAT  ARE  THE  VALUES  FOR  A?  AND  FOR  B“ 

INPUT,  ACON,BCON 

PRINT ,  “WHAT  IS  THE  ANGLE  OF  THE  EARTH  MAGNETIC  FIELD* 

INPUT ,  ANGLE 

BETA  =  ANGLE*PI 

DO  10  V  =  0.01, LIMIT, INC 

AHX  =  CABS ( HX ( FREQ  ,  X ,  Y ,  Z )  > 

WRITE  (1,100)  Y ,  AHX 
WRITE  (If 200 ) 

DO  20  Y  =  0.01, LIMIT , INC 
AHY  =  CAF«S(HY(FREQ,X,Y,Z)  ) 

WRITE  (1,100)  Y,  AHY 
WRITE  (lc 200 ) 

DO  30  Y  =  0»01 , LIMIT , INC 
AHZ  =  CABS ( HZ  <  FREQ  ,  X »  Y ,  Z  )  > 

WRITE  (1,100)  Y ,  AHZ 
WRITE  (1,200) 

DO  40  Y  =  0.01, LIMIT, INC 

AHT  =  CABS(HT (FREQ,X , Y,Z, BETA) ) 

WRITE  (3,100)  Y,  AHT 
WRITE  (3,200) 

DO  50  Y  =  0.01, LIMIT, INC 

WRITE  (3,100)  Y,  HR ( FREO ,X,Y,Z,BETA) 

WRITE  (3,200) 

DO  60  Y  =  0.01, LIMIT, INC 

WRITE  (2,100)  Y , HP ( FREQ , X , Y , Z , BET A ) 

WRITE  (2,200) 

CLOSE  (1) 

CLOSE  (2) 

CLOSE  (3) 

FORMAT  (1PE10.3,1H, ,1PE10.3> 

FORMAT  ( 1 1H1 . E37 , 1 «E37> 

END 
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*  AC  DIPOLE  IN  OU  AS  I  --NEAR  RANGE  *  SUB  TO  AIR 

*  TABLE  3.L3,  KRAI HUMAN 

i: o m p l f: x  f  u N ( : )  h  im  i  i x  <  f  r e o » x » y „ i  > 

COMMON  DEPTH  >■  o  r  tP 
COMPLEX  GAMMA 
R H 0  •-  0 0 K  r  X** 2  i  Y t *  ) 

R  -  SORT  <  RHO -*2  >  2**2) 

GAMMA  t:cjfiRT(l.:MPLX(0.O»FREQ*(i:-I**2)*32.E-7)  ) 

MX  P/  (  2*P 1 1 GAMMA  > 

H  X  ■■  H  X  *  C  t  X  P  (  G  A  M  n  A  *  B  EP  i  H  )  /  <.}<**  3  ) 

MX  -  400  .  >:>  •  1  *HX*3*  (1.0  (  2**2  )  /  (  R**2  )  )  *  (  (  X*Y  >  /  (  RH0**2  )  ) 

R El  UR N 
END 

COMP L E X  P U N C T I  ON  H Y  (  F R E 0  •  X  ,  Y  .  Z  > 

COMMON  DEPTH -P  T  ,  P 

COMPLEX  GArMA 

RHO  -  SORT t X**2  r  Y**2> 

R  --  SOP  TORCH*  *2  4  2**2/ 

GAMMA  -  CSOPT  ( CMP  LX  (0,0,  FREQ*  (  PI**2  >  *32  .  E-7  >  > 

HY  =  P/<2*PI*GAriMA> 

HY  =  400.  *!••■  1  *HY*CEXP  ( 0 AMMA*DEPTH  >  /  <  R**3 ) 

HY  ~  HY  *  •.  2*  <  Y**2  >  -  ( X**2 > - <  3*  < 2**2 ) *  ( Y**2 )  /  < R**2 ))>/( RH0**2 ) 


RETURN 

END 

COMPLEX  FUNCTION  HZ ( FREQ , X » Y , Z ) 

COMMON  DEE ! H , P I rP 

COMPLEX  GAMMA 

RHO  =  SORT ( X**  '  T  Y**2> 

R  =  SORT (RHO** 2  T  2**2) 

GAMMA  =  COURT  ( 1  Mt:'L  X  (  0 . 0  r  FREQ*  ( P  1**2  )  *32 .  E-7  )  ) 

HZ  -  3  «  *  P  /  '  2  *  P  I  *  G  A  M  M  A  *  G  A  h  M  A  ) 

H  7.  «  H  7.  *  Y  *  C  E  X  P  (  G  A  M  M  A  *  D  E  PTH>/<  R  **S) 

HZ  “  400  .  *’• :  YH.  X  <  1  .  O  r  ( GAMMA*Z  > ~5*  (  2**2  )  /  <  E**2  )  ) 

RETURN 

END 

COMPLEX  FUNCTION  H T ( FREQ  r  X  t  Y  ?  Z  r  BETA ) 

COMPLEX  MX.-  Hr,  HZ 
SBETA  -  SIN (BETA) 

CBETA  -  COS (BETA) 

HI  =  HX  ( FREQ  r  X  *  Y  ,  2  >  *CBET A/  •.  .  '  V>  -  0  ♦  823*HZ  ( FREQ  r  X » Y  *  Z ) 

HI  a  HT  t  MY  ( FREQ  rX,Yj7>  *?*£-:  f  A/l .  76 

RETURN 

END 

FUNCTION  HP ( FREQ ?  X »  Y . Z  r BETA ) 

COMPLEX  HX  r  HY -  HZ  r  HT 

HP  =  ATAN(AIMA0  LIT (FREQ fX.YrZr BETA) ) /REAL ( HT ( FREQ ,XfYrZ, BET A ) ) ) 


RETURN 

END 

FUNCT ION  HR ( FREQ  *  X  »  Y  r  Z  r BETA ) 
COMPLEX  HXrHY.HZ 
SBETA  -  SIN (BETA) 

CBETA  -  COS ( BETA ) 
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END 


HR  =  REAL(HX(FREG,X, Y,Z> >#CBETA/1 . 76  -  0.823*REAL(HZ(FREQ,X»Y,Z) ) 

HR  =  HR  +  REAL(HY(FREQ,X,Y,Z) >*SBETA/1.76 

RETURN 

END 

COMPLEX  HX,HY,HZ,HT 
COMMON  DEPTH, P I, P 
DEPTH  =  -76.2 
X=-39 . 5 
PI  =  3.14159 
P  =  50.0 

OPEN  <1, 'AC1FDJ ,ACCESS='ASCII' > 

OPEN  ( 2 , ' AC1PFD' ,ACCESS='ASCII ' > 

OPEN  ( 3f ' AC1TFD *  , ACCESS=  * ASCI I  *  > 

FR£G=1 . 0 

PRINT , 'WHAT  IS  THE  HEIGHT' 

INPUT,  7. 

PRINT,  'WHAT  IS  THE  LIMIT  FOR  Y  AND  ITS  INCREMENT' 

INPUT,  LIMIT , .iX 

PRINT,  'WHAT  IS  THE  ANGLE  OF  EARTH  MAGNETIC  FIELD' 

INPUT,  ANGLE 

BETA  =  ANGLE  *  PI 

DO  10  Y  =  0.01, LIMIT, INC 

AHX  =  CABS ( HX ( FREQ ,X,Y»Z)> 

WRITE! 1 ,100)  Y,  AHX 
WRITE  (1,200) 

DO  20  Y  =  0.01, LIMIT, INC 
AHY  =  CABS  <  HY  <  FREQ ,X»Y»Z)) 

WRITE  (1,100)  Y,  AHY 
WRITE  (1,200) 

DO  30  Y  =  0.01, LIMIT, INC 
AHZ  =  CABS ( HZ ( FREQ ,X,Y,Z)) 

WRITE  (1,100)  Y,  AHZ 
WRITE  (1,200) 

DO  40  Y  =  0.01, LIMIT, INC 

AHT  =  CABS(HT ( FREQ, X,Y,Z, BETA) ) 

WRITE  (3,100)  Y,  AHT 
WRITE  (3,200) 

DO  50  Y  =  0.01, LIMIT, INC 

WRITE  (3,100)  Y,  HR ( FREQ , X , Y , Z , BET A ) 

WRITE  (3,200) 

DO  60  Y  =  0.01, LIMIT, INC 

WRITE  (2,100)  Y,  HP ( FREQ , X , Y , Z , BET A ) 

WRITE  (2,200) 

CLOSE  (1) 

CLOSE  (2) 

CLOSE  (3) 

FORMAT  ( 1PE10 .3,1H,,1PE10.3) 

FORMAT  ( 1 1H1 . E37 , 1 .E37) 
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AC  DIPOLE  IN  QUASI-NEAP  RANGE 
TABLE  3 . 15 f  KRAICKMAN 
COMPLEX  FUNCTION  HX ( FREQ r X r Y , Z ) 

COMMON  DEPT  Hr  PI. P 

COMPLEX  GAMMA 

RHO  =  SGRT(X**2  +  Y*YL; 

GAMMA  =  CSGRHCMPLX(0.0fFREQ*<PI**2>*32.E-7>  > 

HX  *  <3*P)/(2*PX*GAMMA> 

HX  =  HX*CEXP(GAMMA*DEPTH>/(RH0**3> 

HX  -  400  ♦  *P I ¥HX#X* Y /  ( RH0#*2 ) 

RETURN 

END 

COMPLEX  FUNCTION  HY ( FREQ r Xr Y , Z ) 

COMMON  DEPTH t PI rP 

COMPLEX  GAMMA 

RHO  a  SGR!(X**2  +  Y**2) 

GAMMA  =  CSQRT(CMPLX<0.0»FREQ*(PI**2>*32.E-7> ) 

HY  -  !■'/  (  7Y.pl  2GAMMA) 

HY  =--  400.  *F  I*HY*CEXP<GAMMA*BEPTH)/<RH0**3) 

HY  =:  HY  *  (  <2*(  Y**2>  >/<RH0#*2'  •  ( X**2 > / ( RH0**2  )  > 

RETURN 

END 

COMPLEX  FUNCTION  HZ < FREQ , X t Y , Z ) 

COMMON  DEPTH r PI fP 

COMPLEX  GAMMA 

RHO  =  SGRT(X*#2  +  Y**2) 

GAMMA  a  CSQRT(CMPLX<0,0rFREQ*<PI**2>*32.E-7)> 

HZ  =  3.*P/(2*PI*GAMMA*GAMMA> 

HZ  =  HZ*Y*CEXP<GAMMA*DEPTH>/<RH0**5> 

HZ  *  400.*PI*HZ*<1.0  +  GAMMAYZ) 

RETURN 

END 

COMPLEX  Fi  iv-r  T  TON  HI  ( FREQ  r  X  r  Y .  Z  *  BET  A ) 

COMPLEX  Hx - HY r H7 
SBETA  =  SIN  (BETA) 

CBETA  =  COS (BETA) 

HI  =  HX(FRE0rX.Y*Z)*CBETA/1.7A  -  0 . S23*HZ < FREf  . X , Y f Z ) 

HT  =  HT  +  HY(FREQ»XfYrZ)*SBETA'1.76 

RETURN 

END 

FUNCTION  HP  ( FREQ  r  X  r  Y » Z  •  BE  •'£. ) 

COMPLEX  HX « HY  *  HZ  ,*  H  T 

HP  —  AT  AN  ( A I  MAG  ( H  7  ( FREQ » X  r  Y  r  Z  r  BET  A  )  )  /REAL  ( HT  ( FREQ  rXrY  fZj  BETA ) )  ) 

RETURN 

END 

FUNC  T  ION  HR  <  -  RE \Q .  X  r  Y  r  Z  r  BETA ) 

COMPLEX  y f HZ 

SBETA  -  SIN(BETA) 

CBETA  =  COS (BETA) 

HR  =  REV_(HX(FRE0rX»Y*Z))*CBETA/1.76  -  0 . 823*REAL ( HZ (FREQ * X » Y , Z 
HR  =  HR  i  REAL ( HY ( FREQ r X •  Y  ?  Z  ) )YSBETA/1.76 

RETURN 


0»S23*HZ(FREf  .X»Y»Z) 
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END 

COMPLEX  HX*HY*HZ*HT 
COMMON  DEPTH * PI *P 
DEPTH  *  -76*2 
X  =  152.4 
PI  =  3.14159 
P  =  50.0 

OPEN  < 1 * 'ACP2FD* * ACC£SS=* ASCI I •  ) 

OPEN  <  2 » • ACP2PFD * *  ACCESS= ‘ASCII*) 

OPEN  <  3 * " ACP2TFD ‘ *  ACCESS^ • ASC 1 1 ■ ) 

FREQ  =  1. 

PRINT * ‘WHAT  IS  THE  HEIGHT* 

INPUT*  Z 

PRINT *  *WHAT  IS  THE  LIMIT  FOR  Y  AND  ITS  INCREMENT* 
INPUT *  LIMIT * INC 

PRINT *  *WHAT  IS  T.HE  ANGLE  OF  THE  EARTH  MAGNETIC  FIELD* 

INPUT*  ANGLE 

BETA  =  ANGL£*PI 

DO  10  Y  =  0.01* LIMIT* INC 

AHX  a  CABS < HX ( FREQ *X»Y»Z>> 

10  WRITE! 1*100)  Y*  AHX 
WRITE  <1*200) 

DO  20  Y  =  0.01*LIMIT*INC 
AHY  =  CABS  <  HY  t  FREQ  *X*Y*Z)) 

20  WRITE  <1*100)  Y,  AHY 
WRITE  <1*200) 

DO  30  Y  =  0.01 *LIMIT *INC 
AHZ  =  CABS<HZ<FREQ*X* Y  *Z) ) 

30  WRITE  <1*100)  Y*  AHZ 
WRITE  <1*200) 

DO  40  Y  =  0,01 *LIMIT » INC 

AHT  =  CABS  <  HT  <  FREQ  »X*Y*Z»BETA)) 

40  WRITE  <3*100)  Y*  AHT 
WRITE  <3*200) 

DO  50  Y  =  0,01*LIMIT*INC 
50  WRITE  <3*100)  Y*  HR<FREQ*X* Y*Z*BETA) 

WRITE  <3*200) 

DO  60  Y  =  0.01*LIMIT»INC 
60  WRITE  <2*100)  Y*  HP < FREQ *X*Y*Z*BETA> 

WRITE  <2*200) 

CLOSE  <1) 

CLOSE  <2) 

CLOSE  <3) 

100  FORMAT  < 1PE10.3* 1H*  * 1PE10.3) 

200  FORMAT  < 1 1H1 ,E37* 1 . E37) 

END 
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*  POINT  AC  SUB  TO  SUB 

*  TABLE  3.16 

COMPLEX  FUNCTION  HX < FREQ r X r Y , Z ) 

COMMON  DEFTHrPIrP 

COMPLEX  GAMMA1 

RZERO  =  SGRT(X**2  I  Y**2) 

GAMMA 1  =  CSGRT(CMPLX<0.0,FREQ*<PI**2)*32.E~7) ) 

HX  =  P/(2*PI#GAMMAi> 

HX  ==  HX*CEXP(GAMMA1*<  DEPTH  F  Z )  )  /  ( RZER0**3 ) 

HX  =  400,*PI*HX*3*(  (X*Y)/(RZER0**2)  ) 

RETURN 

END 

COMPLEX  FUNCTION  HV ( FREQ » X » Y r Z ) 

COMMON  DEPTH  f  PI  fP 

COMPLEX  GAMMA1 

RZERO  =  SQRT(X**2  +  Y**2) 

GAMMA 1  =  C5QRT(CMPLX<0.0,FREQ*<PI**2)*32.E-7> ) 

HY  =  P/(2*PI*GAMMA1> 

HY  =  400 ♦ *PI #HY#CEXP ( GAMMA1# ( DEPTH  -F  Z ) ) / ( RZER0**3 ) 

HY  =  HY  *  <2*<Y**2>  -  <X**2>  >/<RZER0**2> 

RETURN 

END 

COMPLEX  FUNCTION  HZ < FREQ  *  X ? Y *Z > 

COMMON  DEPTH r PI »P 

COMPLEX  GAMMA 1 

RZERO  =  SQRT<X2*2  +  Y**2) 

GAMMA 1  =  CSQRT(CMPLX(0.0,FREQ*(PI**2>*32.E-7) ) 

HZ  =  3.*P/(2*PI*GAMMA1*GAMMA1) 

HZ  =  HZ*Y*CEXP( GAMMA!* (DEPTH  +  Z ) ) / ( RZER0**5 ) 

HZ  =  400.*PI*HZ 

RETURN 

END 

COMPLEX  FUNCTION  HT ( FREQ » X r Y , Z » BETA ) 

COMPLEX  HX  r HY » HZ 
SBETA  =  SIN<  BEI A ) 

C  BET  A  =  COS (BETA) 

H  r  =  HX  ( FREQ  rXfYfZ)  &CBETA/1 . 76  -  0 . 823*H7.  ( FREn»  X*  Y  *  Z ) 

HI  =  HT  +  HY ( FREQ » X  r Y  r Z  >#SB£TA/1 . 76 

RE  TURN 

END 

FUNCTION  HP ( FREQ r X » Y * Z r BE^A ) 

COMPLEX  HX y  HY  r  HZ  f  HT 

HP  =  ATAN ( A I MAG ( HT <  FREQ  ?  X  j  Y .  Z  r  BET  A )  > /REAL  <  HT ( FREQ r  X » Y  r  Z?  BETA >  > ) 

RETURN 

END 

FUNCTION  HR ( FREQ rXrYrZ, BETA ) 

COMPLEX  KXrHY'KZ 
SBETA  =  SIN (BETA) 

C3ETA  =  COS ( BETA ) 

HR  =  REAL ( HX < FREQ* X r Y » Z ) )#CBETA/1 . ~*6  -  0 . 823*R£AL < HZ< FREQ t X r Y «■ Z ) 

HR  *  HR  +  REAL < HY( FREQ, X f Yr Z ) )*SBETA/ 1.76 

RETURN' 
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30 


40 


50 


60 


100 

200 


END 

COMPLEX  HX*HY*HZ»HT 
DEPTH  =  -76.2 
COMMON  DEPTH * PI *P 
X  =  152.4 
PI  =  3.14159 
P  =  50.0 

OPEN  < 1 *  * ACP2SFD *  *  ACCESS=  *  ASC II*) 

OPEN  (  2  *  " ACP2SPFD *  * ACCESS=  "ASCII* > 

OPEN  <  3  *  * ACP2STFD *  *  ACCESS= *  ASC 1 1  *  > 

FREQ  =  1. 

PRINT * *WHAT  IS  THE  HEIGHT* 

INPUT *  Z 

PRINT *  'WHAT  IS  THE  UPPER  LIMIT  FOR  Y  AND  ITS  INCREMENT* 
INPUT *  LIMIT * INC 

PRINT * *WHAT  IS  THE  ANGLE  FOR  EARTH  MAGNETIC  FIELD* 

INPUT*  ANGLE 

BETA  =  ANGLE*PI 

DO  10  Y  =  0«01*LIMIT»INC 

AHX  =  CABS (  HX  ( FREQ  »X*Y*Z)) 

WRITE< 1 * 100)  Y *  AHX 
WRITE  ( 1 *200) 

DO  20  Y  =  0*01 » LIMIT t INC 
AHY  =  CABS  <  HY  <  FREQ  *X»Y»Z>) 

WRITE  (1»100)  Y»  AHY 
WRITE  < 1 »200) 

DO  30  Y  =  0 . 01 » LIMIT t INC 
AHZ  =  CABS  <  HZ  <  FREQ  »X»Y»Z)) 

WRITE  <lflOO)  Y,  AHZ 
WRITE  < 1 »200) 

DO  40  Y  =  0.01 » LIMIT » INC 

AHT  =  CABS(HT(FREQ»X»Y»ZfBETA)> 

WRITE  <3r 100)  Y»  AHT 
WRITE  <3»200) 

DO  50  Y  =  0.01 rLIMIT *  INC 

WRITE  (3*100)  Y f  HR(FREQ* XrYrZ»BETA) 

WRITE  (3*200) 

DO  60  Y  =  0.01  * LIMIT* INC 

WRITE  (2*100)  Y*  HP ( FREQ *X*Y*Z»BETA) 

WRITE  (2*200) 

CLOSE  (1) 

CLOSE  (2) 

CLOSE  (3) 

FORMAT  ( 1PE10. 3* 1H*  *1 PE 10. 3) 

FORMAT  ( 11H1 .E37* 1 .E37) 

END 
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it 

;  *  INFINITFSIMAl  DC  HE  D 

»• 

-8 

1  FUNC  T I  ON  HX  *  X  *  Y  *  7  '■ 

!  i  OMrlUN  i.-F P  I  H  *  F 

a 

|  V  7  FIFTH 

ji 

1  KHO  -  sNklvX**:’  }■  Ts*7' 

i 

1  K’lJNf  -  *,URl  <RH0#*2  i  •  H#*2 ) 

1 

B  HX  -  t  F  *X>rT  '•  • : iKHO*#!’  • 

n 

|  Rf ’ URN 

1  E  NU 

n 

p 

■;  F  UNL  s  I  ON  HY(X*Yfi  ' 

S 

1  common  hepth*p 

-j! 

1  ZH  =  Z  -  DEPTH 

|  PHD  -  SORT  (X**2  +  Y**2> 

H 

1  RONE  =  SQF7(RHG**2  +  ZH**2> 

a 

7 


I  HY  =  J  OO  ,  *P*HY  '*  (  FH0*#2  ) 

S 

I  rv f- 1  s jri'f- 

|  Fr.-n 

|| 

|  F  UN f'  1 1  ON  HZ  ( X  r  Y  *  Z ) 

§  COMMON  D.f  f  -  rH*F- 

-7» 

fl 

I  Zh  -  Z  ~  EH 

Tjl 

1!  RHO  -  SURUX**2  +  Y**2> 

RONE  -  SORT'.KHO**?  +  ZH**2> 

|  HZ  -  100*P*Y  RONE*  <  RONE**?  >  > 

|  RETURN 

5 

r  END 

1 

|  FUNCT  ION  HI ( X »  Y • Z  r  BETA) 

f  SBETA  =  PIN < BETA) 

I  CFiEl  A  =  COS  <  BETA ' 

1 

p:  RETURN 

1 

I;  E.ND 

H  COMMON  DL.  :(!»P 

1 

1 

1  Dr'FT  H=-76 . 2 

j|  •  =-152.4 

B  p  =  5o.O 

E’ 

1  PI  =  3.14159 

I  OPEN  fit’  F»C1  FIELD*  *  ACCFSS=  *  ASC7  I  °  ) 

'  'Ej 

ri 

i  ENDULE  ' 

,?■ 

I  PRINT  *  *  WHAT  IS  THE  HF  T Gu 1 

1  INPUT*  7 

•  V 

|  INPUT  *  i  I H  T  f  *  i  fH‘ 

a 

|  PRINT?  ’CHAT  I!  THE  ANGLE  OF  THE  EARTH 

MAGNETIC  FIELD*  § 

I  INPUT*  AN*'  - 

1 

W  RE  •  A  -  A*-*  *!F*P1 

j| 

1  I*r«  10  --  -L IMIT  ?!.IMIT  y  INC 

I  10  WRITF  (1?100)  Y  *  HX  ( X  ?  Y  *  Z  .) 

1§ 

1  WRITE  ■ 1*400? 

1  DC  20  Y  =  -LIMIT. LIMIT* INC 

1  20  WPI'-  ■ 1 ? 100)  Y*  HY(X.YfZ) 

a 

|  07/23/00 

ii:09:si  |j 

1  2  AI6 

I 

Hl,l  -  . _ - 

 ...  -  J 

DCP 


(PaSe  2) 


WRITE  (If 400) 

DO  30  Y  =  -LIMITfLIMITf INC 
30  WRITE  (lflOO)  Y»HZ(XfYfZ> 
WRITE  (If 400) 

DO  40  Y  =  -LIMITfLIMITf INC 
40  WRITE  (lflOO)  Yf HT(Xf Yf Zf BETA) 
WRITE  (If 400) 

CLOSE  (1) 

)0  FORMAT  < 1PE10. 3f lHf  1 1PE10. 3 ) 

)0  FORMAT  (HHl.E37fl.E37) 

END 
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*  MODIFIED  IMAGE  FOR  FINITE  LENGTH  DC  HED 
FUNCTION  HX(XrYfZ) 

COMMON  BCON r DEPTH  *  CURENT  f  LENGTH 

REAL  LENGTH 

REAL  K21 r K22 

XLP  =  X  -f  LENGTH/2 

XLM  =  X  -  LENGTH/2 

ZB  =  Z  -BCON*DEPTH 

Xl.PY  =  XLP**2  +  Y**2 

XLMY  =  XLM**2  +  Y**2 

ZBY  -  ZB**2  +  Y**2 

SORT ( Xl„P**2  +  ZBY) 

SGRT(XLM**2  +  ZBY) 

(1  -  ZB/K22  > /XLMY 
<1  -  ZB/K21 )/XLPY 
100*Y*CURENT*(HX1  -  HX2) 


K21  = 
K22  « 
HX1  = 
HX2  = 
HX  = 


RETURN 

END 

FUNCTION  HY ( X  f  Y » Z ) 

COMMON  BCON  r  DEPTH  r CURENT  >  LENGTH 
REAL  K21 » K22 
REAL  LENGTH 
XLP  =  X  +  LENGTH/2 
XLM  a  X  -  LENGTH/2 
ZB  =  Z  -BCON*DEPTH 
XLPY  =  XLP**2  +  Y**2 
XLMY  =  XLM**2  +  Y**2 
ZBY  a  ZB**2  +  Y**2 
K21  a  SGRT(XLP**2  +  ZBY) 
SQRT(XLM**2  +  ZBY) 

( XL.P/K21  -  XLM/K22 >*ZB/ZBY 
XLP*  < 1  -  ZB/K21 ) /XLPY 
XLM* ( 1  -  ZB/K22 ) /XLMY 
1 00 ♦ *CURENT*  < -HY 1  +  HY2 


a 


K22  = 
HY1  = 
HY2  = 
HY3  = 
HY  = 


HY3 ) 


RETURN 

END 

FUNCTION  HZ ( X  r  Y  » Z ) 

COMMON  BCON t DEPTH r CURENT r LENGTH 
REAL  K21 r N'22 
REAL  LENGTH 
XLP  =  X  +  LENGTH/2 
XLM  a  X  -  LENGTH/2 
ZB  =  Z  -BCON*DEPTH 
ZBY  =  ZB**2  +  Y**2 
K21  =  SORT ( XLP**2  + 

K22  =  SORT (XLM**2  + 

HZ  =  1(>0*Y*CURENT*(  XLP/K21 
RETURN 
END 

FUNCTION  HT(X»YfZrBETA) 

SBETA  =  GIN (BETA) 

CBETA  =  COS (BETA) 


ZBY) 

ZBY) 


XLM/K22 ) /ZBY 
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HT  =  HX(X»Y*Z)*CBETA/1.76  +  HY(X*Y*Z)*SBETA/1 .76  -  0.823*HZ(X*Y*Z) 

RETURN 

END 

REAL  LENGTH 

COMMON  BCON* DEPTH *CURENT *LENGTH 

DEPTH=-76.2 

CURENT  =1.00 

BCON  =1.00 

LENGTH  =50.0 

X=304.8 

PI  s  3.14159 

OPEN  (1»  * DCF I ELD* * ACCESS="ASCII " ) 

ENDPILE  1 

PRINT *  ‘WHAT  IS  THE  HEIGHT  IN  METERS" 

INPUT *  Z 

PRINT *  "WHAT  IS  THE  LIMIT  FOR  Y  AND  ITS  INCREMENT" 

INPUT *  LIMIT * INC 

PRINT *  "WHAT  IS  THE  ANGLE  FOR  EARTH  MAGNETIC  FIELD" 

INPUT*  A 
BETA  =  A*PI 

DO  10  Y  =  “LIMIT * LIMIT* INC 
10  WRITE  (1*100)  r'?HX<X*Y»Z> 

WRITE  (1*200) 

DO  20  Y  =  ~LIMIT*LIM!T»INC 
20  WRITE  (1*100)  Y*HY(X»Y»Z) 

WRITE  (1*200) 

DO  30  Y  =  -LIMIT*LIMIT*INC 
30  WRITE  (1*100)  Y,  HZ(X*Y»Z) 

WRITE  (1*200) 

DO  40  Y  =  -LIMIT»LIMIT»INC 
40  WRITE  (1*100)  Y*HT(X*Y*Z*  BET  A ) 

WRITE  (1*200) 

CLOSE  (1) 

100  FORMAT  (1PE10.3*1H**1PE10.3) 

200  FORMAT  ( 11H1 .E37* 1 .E37) 

END 
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*  INFINITESIMAL  DC  HED  FROM  SUBSURFACE  TO  SUBSURFACE 

*  FORMULAE  (3.18)  -  (3.20) 

FUNCTION  HX(XrYi’Z) 

COMMON  DEPTH fP 

ZHP  ==  Z  +  DEPTH 

RHO  =  SQRT(X**2  +  Y**2> 

RTWO  =  SQRT(RH0**2  ZHP**2) 

HX  =  (P*X*Y)/(RH0**2) 

HX  =  100.*HX*(ZHP/(RTW0**3>  +  2*(ZHP/RTW0  +  1>/(RH0**2)) 

RETURN 

END 

FUNCTION  HY(X»Y»Z) 

COMMON  DEPTH fP 

ZH  =  Z  -  DEPTH 

ZHP  =  Z  +  DEPTH 

?H0  =  SORT ( X##2  +  Y**2) 

RONE  =  SQRT(RH0**2  +  ZH**2) 

RTWO  »  SGRT(RH0**2  +  ZHP**2) 

HY1  =  ZH/(R0NE**3)  +  ZHP*<X**2>/( (RH0**2)*(RTW0**3) ) 

HY2  =  (ZHP/RTWO  +  1)*(X**2  -  Y**2)/(RH0**4> 

HY  ==  -100,*P*(HY1+HY2> 

RETURN 

END 

FUNCTION  HZ(XfYrZ) 

COMMON  DEPTH rP 
ZH  =  Z  -  DEPTH 
RHO  =  SGRT(X**2  +  Y**2> 

RONE  =  SQRT(RH0**2  +  ZH**2) 

HZ  =  10O*P*Y/<RONE*(RONE#*2>> 

RETURN 

END 

FUNCTION  HT(XrYrZrBETA) 

SBETA  =  SIN (BETA) 

CBETA  =  COS (BETA) 

HT  =  HX(XfYrZ) *CBETA/1 . 76  +  HY ( X r Y r Z ) *SBETA/1 . 76  ~  0 . 823*HZ ( X » Y t Z ) 

RETURN 

END 

COMMON  DEPTH fP 
DEPTH=-76 « 2 
X=1 

P  =  50.0 
PI  =  3.14159 

OPEN  ( 1  * ’ DCSUBFD *  t ACCESS= ‘ ASC II") 

ENDFILE  1 

PRINT r * WHAT  IS  THE  HEIGHT* 

INPUT  i  Z 

PRINT t  -WHAT  IS  THE  LIMIT  FOR  Y  AND  ITS  INCREMENT* 

INPUT »  LIMITtINC 

PRINT.  "WHAT  IS  THE  ANGLE  OF  THE  EARTH  MAGNETIC  FIELD* 

INPUT t  ANGLE 

BETA  a  ANGLE#PI 

DO  10  Ya  -L IMITr LIMIT » INC 
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10 


20 


30 

40 


100 

400 


WRITE  ( 1 , 100 >  Y,HX(X,Y,Z> 

WRITE  ( 1  * 400 ) 

DO  20  Y=  -LIMIT, LIMIT, INC 
WRITE  (1,100)  Y,  HY(X,Y,Z) 
WRITE  ( 1 ,  400  > 

DO  30  Y=  -LIMIT, LIMIT, INC 
WRITE  ( 1 , 100 )  Y,HZ(X,Y,Z) 

WRITE  (1,400) 

DO  40  Y=  -LIMIT, LIMIT, INC 
WRITE  (1,100)  Y ,HT< X, Y ,Z,BETA) 
WRITE  (1,400) 

CLOSE  (1) 

FORMAT  (1PE10»3,1H,,1PE10»3) 
FORMAT  (11H1.E37,1»E37> 

END 
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